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Abstract. 

Some new results on magnetic dipole (Ml) transitions in heavy quarkonium from nonrelativistic 
effective field theories of QCD are briefly reported. This model-independent approach not only 
facilitates a systematic and lucid way to investigate the relativistic corrections, it also clarifies 
some inconsistent treatment in previous potential model approach. The impact of our formalism 
on J/y — > rjc-Y, T(T') — > r\ b y and h c — > XcoY are discussed. 
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Radiative transitions in heavy quarkonium are of considerable experimental and theo- 
retical interest [ 1 ] . On the theory side, it provides us with further insight on the dynamics 
of quarkonium in addition to the knowledge we have gleaned from the spectra of cc and 
bb families. 

^ \ Being an old subject, radiative transitions have been extensively studied in phe- 

nomenological models, notably the potential model approach [I2j]. It is certainly desir- 
able to study them from a model-independent perspective. In this talk I will report such 
P-i! a study based on the effective-field-theory (EFT) approach [3]. Since Ml transition is 

theoretically cleaner and more interesting than El transition, I will focus on the former 
case, despite the fact that the latter is observed more copiously in nature. 

In the nonrelativistic limit, the Ml transition rate between two 5-wave onia takes a 
particularly simple form: 
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where eg and Kq are the electric charge and anomalous magnetic moment of the heavy 
quark, and R n i(r) is the radial Schrodinger wave functions. Since the leading dipole 
operator CFeeg/2mc ■ B em (where cp = 1 + Kq) only flips the spin and doesn't act on 
the spatial degrees of freedom, orthogonality of radial wave functions guarantees that: in 
the allowed transition (n = n'), the overlap integral equals 1; in the hindered transition 
(n 7^ n'), the overlap integral vanishes. 

It is well known that (U]) overpredicts the observed J / y — > r) c y transition rate by a fac- 
tor of 2 ~ 3, with a normal input of m c and K c from perturbative one-loop matching. This 
clearly indicates large relativistic corrections to ©, as is usually confronted in charmo- 
nium system. It has also been speculated that low-energy fluctuations may generate a 
large negative K c , so the discrepancy can be reduced. 



The &(y 2 ) corrections to has been available for a long time from potential model 
approach JEH: 
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where 27„ + 1 counts number of polarizations of the parent onium, and 
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where Vb stands for the static potential, the "+/-" sign in I4 is associated with 3 5i — > ^0 
and 5o — >■ 3 Si, respectively. Notice that I4 accounts for the first-order correction to the 
wave function due to spin- spin potential (other higher dimensional local potentials cease 
to contribute in S-wave transition), thus is only present in hindered transition. 

/5 is a prediction specific to the popular assumption in potential models, where one 
usually decomposes the confining potential into a Lorentz scalar and a vector part, 
Vconf = TfVs + (1 — Tj)Vv- This term constitutes the major uncertainty in the potential 
model predictions, where contradictory claims often appear in the literature |0]. 

Presence of a hierarchy of scales in quarkonium, m, mv, mv 2 , Aqcd, makes the EFT 
approach an ideal tool to analyze this transition process. One first descends from QCD 
to NRQCD by integrating out hard modes (~ m ) ||4|], then descends from NRQCD 
to potential NRQCD (pNRQCD) by further integrating out soft modes (~ mv ) o. 
As a result, the inter-quark potentials appear as Wilson coefficients in pNRQCD, and 
the only dynamical degrees of freedom of pNRQCD are ultrasoft modes (~ mv 2 ) (for 
convenience, one can also incorporate the radiated ultrasoft photon into pNRQCD). 

A primary task of the EFT approach is to validate/invalidate ©• The real strength of 
the EFT approach is, however, that it can further answer the following questions that are 
beyond the scope of potential models: 

• Is it possible that a large correction to Kq due to soft modes arises when one 
descends from NRQCD to pNRQCD? 

• Is it possible to reproduce 1$ in pNRQCD? If so, how to interpret it? 

• Potential model focus exclusively on the QQ Fock-state (with an exception of 
coupled-channel effects which may not be relevant as long as one excludes those 
states close to the open-flavor threshold). pNRQCD allows one to include ultrasoft 



gluons as dynamical degrees of freedom. Therefore, one naturally asks for the 
possibility of large nonperturbative contribution arising from higher-Fock states 
\QQg). This effect is usually referred to as color-octet effect. 

In fact, the answers to these three questions are all negative 0], on which we will 
elaborate shortly. It turns out that the pNRQCD formalism is able to justify © except 
I5. In the so called weak-coupling regime (when mv 3> AqcdX this formula is complete; 
in the so-called strong-coupling regime (when mv ~ AqcdX however, © is incomplete 
and further terms are needed. 

Before going on to explanations, it is useful to first sketch the derivation of 0(v 2 ) 
corrections in the framework of pNRQCD. One obvious source is from the contribution 
of higher dimensional Ml operators. These operators can be identified most conveniently 
by promoting the color gauge group of NRQCD Lagrangian to a larger gauge group 
SU C (3) x C/ em (l). Explicitly, the relevant magnetic operators up to ^(1/m 3 ) read |6]: 
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where = + igT a A a ^ + ieeQA s ™. Various Wilson coefficients can be computed 
through perturbative matching at the hard scale. For example, at one loop accuracy, 
Kq = Cf a s /2n, is about a few percent for charm and bottom. It is well known that 
some of the Wilson coefficients are related with each other because of reparameterization 
invariance (RPI), or essentially Poincare invariance, notably cs = 2c f — 1, cw2 = CW\ — 1 
and c p i p = cp — 1 [6]. When inherited into pNRQCD, one make replacement V — > z'p in 
these operators, and these RPI relations can be utilized to condense the expressions. 

There are other sources of 6{y 2 ) corrections, for instance, multipole expansion of the 
photon field. One interesting contribution, first pointed out by Grotch and Sebastian 
is the Lorentz boost effect due to the final-state recoil. Since the wave function of a 
moving 5-wave state has a non-vanishing overlap with a P-wave, spin-flipped state at 
rest, the Ml transition can be effectively realized by a "El " transition from the parent to 
this small component. Some subtlety arises in this effect, namely the recoil correction 
depends on which "El" operator one uses, i.e., 2eeg/mp-A em or ee^r ■ E em , which 
are connected by a redefinition of pNRQCD field. The solution to this problem is as 
following. The matching from NRQCD generates a new Ml operator at €?(l/m 2 ), which 
is intimately related to the non-local (depending on the center-of-mass momentum) spin- 
orbit potential, with a coefficient proportional to V^ M = — V^/Sr (Gromes relation). 
The contribution of this operator also depends on the field redefinition. However, the 
sum of these two corrections is convention-independent, thus comprises a meaningful 
relativistic recoil effect. 

We are now in a position in discussing general matching of Ml operators from 
NRQCD to pNRQCD. At ff{l), those operators in © are trivially inherited to pN- 
RQCD. In general, after integrating out soft modes, new Wilson coefficients will depend 
on inter-quark separation r. 




FIGURE 1. Typical NRQCD diagrams responsible for matching of the dipole operator at ff(l/m ). 
The left diagram represents the vertex correction, and the right one contributes to the wave-function 
renormalization factor Z s (r) due to soft modes, where the cross-cap implies insertion of a unit operator. 



Obviously the correction to cf can be interpreted as a multiplicative &{\/rrr) match- 
ing coefficient of the leading Ml operator. Dimensional considerations require that this 
correction is function of log(r). We emphasize in passing that it is inappropriate to at- 
tribute this new coefficient to the magnetic moment of an individual quark, since it really 
arises from entangled contributions from both quarks. 

The simple answer, no corrections at all, is basically nothing more than the heavy 
quark spin symmetry and that the o ■ B em operator behaves like a unit operator in 
spatial Fock space. These facts are independent of whether the matching is performed 
perturbatively as in weak-coupling regime, or nonperturbatively, as in strong-coupling 
regime. In an arbitrary NRQCD diagram depicting transition, as shown in Fig. Uk), we 
only need to consider the case where the external ultrasoft photon is attached to one 
of the heavy quark lines, (photon attached to internal light quark loops can be neglected 
when summing over electric charges of three light flavors). Soft gluons attached to heavy 
quark lines must be longitudinal to avoid 1/m suppression. Because all the propagators 
and vertexes are spin-independent except the Ml vertex, such a diagram can factorize 
into the leading Ml operator times minus the wave function renormalization constant 
8Z s (r), which can be extracted from the same diagram with the electromagnetic vertex 
replaced by a unit operator insertion, as shown in Fig.[T]b). Therefore, there is no net 
contribution to this matching coefficient from soft modes. 

The matching at &{ljm) is slightly more complicated. In a NRQCD diagram with 
a insertion of the leading dipole opeartor, a transverse gluon is allowed to end on the 
heavy quark line with a spin-independent p • A vertex (ending on a a ■ B vertex doesn't 
contribute). Repeating the previous argument, one doesn't obtain any new Ml operator 
from these sources. However, there are a new class of diagrams, in which electromag- 
netic coupling is embedded in the covariant derivative D, inducing an effective magnetic 
operator at &{\/m). This is just the one accompanying with the nonlocal LS-potential 
as mentioned before. It is important to note that there is no room at this order for the 
operator of the form V^/m 2 o ■ B em to arise. Therefore one is forced to regard /5 in © as 
an artifact of potential model, which cannot bear any physical significance. 

The matching at <^(l/m 2 ) becomes much more cumbersome. We will not dwell on 
the explicit expressions for these case. However, it should be kept in mind that they 
may be neglected in the weak-coupling regime because of suppressions by higher pow- 
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FIGURE 2. pNRQCD diagrams for color-octet contribution to radiative transition. Single and double 
lines correspond to the singlet and octet fields. The singlet-octet-gluon vertex is of chromo-El type. 



ers of a s ~ v. In the strong-coupling regime, since a s , ~ 1, these operators might be as 
important as other @{y 2 ) corrections. These corrections involve some unknown nonper- 
turbative Wilson loop amplitudes, so that predictive power is unavoidably damaged. 

The color-octet effect to radiative transitions was first envisaged by Voloshin long 
time ago, but without a detailed study [7]. pNRQCD allows a systematic treatment of 
this effect. The corresponding diagrams are shown in Fig. El It is not surprising that the 
color-octet effect is nearly absent in Ml transition, because of the exactly same reason 
as before. It is interesting to note this bears some resemblance as absence of leading 
nonperturbative correction to B — > D* form factor at the zero recoil, known as Luke's 
theorem [8]. Here we only give a heuristic argument based on the quantum-mechanical 
perturbation theory. Treating the chromo-El interaction as perturbation up to second 
order, one can schematically express a "true" quarkonium state as superpositions of the 
following Fock components (stripping off spin d.o.g.): 

| AO = fifii (/»)) + 1 QQ % g) + £ | QQi (m) ) • • • , (5) 

where \QQ\in)) is the unperturbed state, and Z" s is the wave function renormalization 
factor for state n due to ultrasoft gluons. Spin-independence of chromo-El interaction, 
plus Ml operator being a unit operator in coordinate space, imply that the full Ml tran- 
sition amplitude is nothing but the inner product (N'\N). Since orthogonality condition 
is preserved in perturbed states, it is equal to (n'\n), the color-octet contribution thus 
vanishes. We emphasize that sizable color-octet effect might arise in El transition, since 
the El operator is no longer a unit operator in coordinate Fock space. 

Finally we turn to the phenomenological implication of ©, with the understanding 
that we restrict only to the weakly-coupled system for consistency and the unphysical I5 
term has been dropped. Thus far, the only observed Ml transitions are 7/ y(vO TlcY, 
and upper bounds on T'(T") — > r\bj have recently been set by CLEO [10]. Empirically, 
T(15) and rjt, are believed to lie in the weak-coupling regime, whereas J/\jf, T\ c fit 
in this regime to a less extent. T' and T" are usually regarded as strongly-coupled 
system. As for y/ and T]' c , they are too close to the open-flavor threshold and cannot 
be correctly described by current formulation of pNRQCD, therefore we exclude them 
in our analysis. 

For a weakly-coupled system, the dynamics is largely governed by the perturbative 
static potential, i.e., Vq ~ —CfOc(ii) /r, where the natural choice of /1 is around the 
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FIGURE 3. r[//y/ — > 77 c y] (in keV) vs. ;ii (in GeV). Dashed line is the prediction in NR limit, solid 
line includes Giy 1 ) corrections, and the band represents the measured width with error taken from |9]. 



typical three-momentum scale. One also needs to specify the quark mass in ©. A naive 
input is to use m, half of the center-of-gravity ground state mass. However, this simple 
procedure may induce an error of order v 2 , which shouldn't be neglected according to 
power counting. A more consistent way is to choose the 15 mass m, which is defined 
implicitly through m = m- (10 \p 2 /2m\ 10) « m(l -C 2 7 a. s ( i u) 2 /8). 

Many of terms in © are practically negligible. Since Kq retains its NRQCD value and 
is only a few per cent for charm and bottom, we may simply put it to be zero in those 
&(y 2 ) terms (as a result, It, can be dropped in both allowed and hindered transitions). In 
fact, only I2 in accounts for genuine &{v 2 ) correction for allowed transition, and we 
end up with a simple expression 

, „ 2C 2 F a 2 (uY 
\+2k c . (6) 

Fig. |3] shows a comparison between this formula and the data. It seems that © is 
compatible within the error to the data when /! is about 0.8 GeV. Note this value is 
consistent with the empirical mv value for J /yr. Therefore, our reasonable success in 
describing J/\\f — » r\ c y may be viewed as a posteriori support for the weak-coupling 
assignment. However, rather sharp ji dependence may suggest that J/\]/ is not far from 
the strong-coupling regime. One can employ © to T — > r\},y with more confidence, 
and finds a smaller 6{y 2 ) correction with flatter scale dependence. Unfortunately, very 
narrow width of 2 ~ 3 eV makes this transition unlikely to be detected. 

For hindered transitions, e.g. T'(T") — > rjhY, the experimental upper bounds are al- 
ready rather tight and many model predictions have been ruled out [ 10]. To be conserva- 
tive, one would not expect © to reliably describe these transitions if excited bottomonia 
are indeed in the strong-coupling regime. Nevertheless, it may not be too optimistic to 
treat, at least T', as being in weak-coupling regime. Proceeding along this line, we find 
that the I4 term is dominating over I\ and I2, and the latter two nearly cancel each other. 
As a result, the predicted width is about an order-of-magnitude larger than the experi- 
mental upper bound! Even though there are lots of uncertainty associated with I4, this 
alarmingly large discrepancy seems to indicate that the weak-coupling assignment of T' 
is problematic and a strong-coupling analysis might be more appropriate. 



We end with a brief discussion on the P-wave Ml transitions, which have received 
less attention in the past, due to lack of phenomenological impetus. Recent discovery 
of the h c state 111 ill may arouse the interest to study these processes. Color-octet effect 
again vanishes because of the exactly same reason. It is straightforward to apply the 
pNRQCD formalism to derive the 6{y 2 ) corrections. For simplicity, we just quote the 
weak-coupling formula for the allowed transition: 
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where cj = 1/2, 1,4/5 for 7 = 0, 1,2. For the n P\ — > n ^Pjj transition, one simply mul- 
tiplies the right side of © by a statistical factor of (27+ l)/3. It should be understood 
that this formula may be of limited use, since P-wave onia may necessarily live in the 
strong-coupling regime. 

The fine splittings between h c and XcO> Xci an d Xa are about 110, —13, —32 MeV, 
respectively. It seems that only h c — > XcoY nas a serious chance to be observed, with a 
width comparable to T[J/yf — ► r\ c y]. Future experimental input will enable us to infer 
the size of relativistic corrections, thus enriching our understanding of P-wave onia. 
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